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ABSTRACT 

We present [Fe/H] abundance results that involve a sample of stars with a wide 
range in luminosity from luminous giants to stars near the turnoff in a globular cluster. 
Our sample of 25 stars in M71 includes 10 giant stars more luminous than the RHB, 
3 horizontal branch stars, 9 giant stars less luminous than the RHB, and 3 stars near 
the turnoff. We analyzed both Fe I and Fe II lines in high dispersion spectra observed 
with HIRES at the W. M. Keck Observatory. We find that the [Fe/H] abundances 
from both Fe I and Fe II lines agree with each other and with earlier determinations. 
Also the [Fe/H] obtained from Fe I and Fe II lines is constant within the rather small 
uncertainties for this group of stars over the full range in T e f f and luminosity, suggesting 
that NLTE effects are negligible in our iron abundance determination. In this globular 
cluster, there is no difference among the mean [Fe/H] of giant stars located at or above 
the RHB, RHB stars, giant stars located below the RHB and stars near the turnoff. 

Subject headings: globular clusters: general — globular clusters: individual (M71) - 
stars: evolution - stars:abundances 



1. INTRODUCTION 

Abundance determinations of stars in Galactic globular clusters can provide valuable informa- 
tion about important astrophysical processes such as stellar evolution, stellar structure, Galactic 
chemical evolution and the formation of the Milky Way. Surface stellar abundances of C, N, O, 
and often Na, Mg, and Al are found to be variable among red giants within a globular cluster. The 
physical process responsible of these star-to-star element variations is still uncertain (see Cohen et 
al. 2001, Paper I). 
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Of particular importance to the present study are the results of King et al. (1998), who found 
that [Fe/H] among M92 subgiants is a factor of two smaller than [Fe/H] from red giants in the 
same cluster (Cohen 1979; Sneden et al. 1991). If this result is not due to systematic differences 
arising from the analysis procedures of the different groups that handle the red giant branch (RGB) 
and the subgiant samples, then the results of King et al. (1998) would suggest some modification 
of photospheric Fe abundances which would be quite difficult to explain. Clearly the present day 
clusters stars are incapable of modifying their own Fe abundances through nucleosynthesis. Yet 
evolution-driven changes in [Fe/H] could be possible if the outer envelopes of the subgiants were 
somehow infused with Fe-poor external material. Thus as evolution progesses up the RGB the 
deeping convective envelope would dilute the photospheric [Fe/H] with more Fe-rich "unmodified" 
material from the interior causing an increase in [Fe/H] with luminosity. We note the unlikelyhood 
of this scenario as Fe, unlike the lighter elements, is not involved in normal mass transfer processes 
that might ocurr in binary systems or in normal stellar winds, and is only produced in supernovae. 
Moreover, it is difficult to understand how cluster stars could be contaminated by material with 
less Fe, presumably at a time after the present subgiants had formed. 

In order to study the origin of the star-to-star abundance variations and to address the issues 
raised by King et al. (1998), we have started a program to determine chemical abundances of the 
nearer galactic globular cluster stars. In this paper, we present our results for the iron abundance 
of M71, the nearest globular cluster reachable from the northern hemisphere. Similar programs are 
underway at ESO by Castilho et al. (2000) and Gratton et al. (2001) taking advantage of the fact 
that the nearest globular cluster accessible from a southern site (NGC 6397) is 1.8 times close than 
M71 and has a lower reddening. Our M71 sample includes stars over a large range in luminosity: 19 
giant stars, 3 horizontal branch stars, and 3 stars near the main sequence turnoff, in order to study 
in a consistent manner red giants, horizontal branch stars, and stars at the main sequence turnoff. 
Details on the star sample, observations, data reduction and determination of stellar parameters 
are described in Paper I. Previous high dispersion abundance analysis for M71 involve studies of 
red giants only, resulting in [Fe/H] of -0.70 (Cohen 1983), -0.80 (Gratton et al. 1986), between 
-0.6 and -1.0 (Leep et al. 1987), and -0.79 (Sneden et al. 1994). 



2. ANALYSIS 

We begin our analysis with Fe, as many Fe lines are identified in our HIRES spectra over a 
wide range of excitation potentials and line strengths, as well as over two ionization states (Fe I and 
Fe II). The iron abundance analysis is done using a current version of the LTE spectral synthesis 
program MOOG (Sneden 1973). A line list specifying the wavelengths, excitation potentials, gf- 
values, damping constants, and equivalent widths for the observed Fe I and Fe II lines is required. 
The provenance of the gf-values and the measurement of equivalent widths are discussed below. The 
damping constants for all Fe I and Fe II lines were set to twice that of the Unsold approximation for 
van der Waals broadening following Holweger et al. (1991). The use of the Blackwell approximation 
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for the damping constants gives the same result within the errors, when comparing [Fe/H] obtained 
with the "good line set" of Fe I (see definition of line sets on Sec. 2.2). 

In addition, a model atmosphere for the effective temperature and surface gravity appropriate 
for each star and a value for the microturbulent velocity are also required. We use the grid of 
model atmospheres from Kurucz (1993b) with a metallicity of [Fe/H] = —0.5 dex, based on earlier 
high dispersion abundance analysis of M71 red giants (Cohen 1983; Gratton et al. 1986; Leep et 
al. 1987; Sneden et al. 1994). The final result for [Fe/H] is not sensitive to small changes in the 
metallicity of the model atmosphere. In particular, the error introduced by using a model with 
[Fe/H] = —0.5 instead of —0.7 dex is very small (see Tables 2 & 3). The effective temperatures and 
surface gravities are derived from the photometry of the stars as described in Paper I. The error in 
the photometric T e ff is 75 K for giants and 150 K for the dwarfs and the error in the photometric 
log(g) is 0.2 dex (Paper I). The microturbulent velocity is derived spectroscopically (see below). 
The stellar parameters are listed in Table 1. 

2.1. Transition Probabilities 

Transition probabilities for the Fe I lines were obtained from several laboratory experiments, 
including studies of Fe I absorption lines produced by iron vapor in a carbon tube furnace (Blackwell 
et al. 1979, 1982a,b, 1986) (Oxford Group), measurement of radiative lifetimes of Fe I transitions by 
laser induced fluorescence (O'Brian et al. 1991; Bard et al. 1991; Bard & Kock 1994), Fe I emission 
line spectroscopy from a low current arc (May et al. 1974), and emission lines of Fe I from a shock 
tube (Wolnik et al. 1971). We also considered solar gf-values from Thevenin (1989, 1990) when 
needed. 

We compare the gf-values obtained by the different experiments in an attempt to place them 
onto a common scale with respect to the results from O'Brian et al. (1991), who provided the 
longest list of gf-values. We considered for the comparison only the set of lines present in our 
data, which have a wavelength coverage roughly from 5380 A to 7900 A. We found that the values 
of O'Brian et al. (1991) and of the Oxford Group were on the same scale; the mean difference in 
log(gf) between the two experiments is 0.02 ±0.01 for 21 lines in common. Similar results are found 
when considering 21 lines in common between O'Brian et al. (1991) and Bard et al. (1991) (see 
also Bard & Kock 1994), where the mean difference in log(gf) is 0.02 ±0.04. Considering 34 lines 
in common between O'Brian et al. (1991) and May et al. (1974), the mean difference in log(gf) is 
0.03 ± 0.03. When comparing the 14 lines in common between O'Brian et al. (1991) and Wolnik 
et al. (1971) we found a mean difference in log(gf) of —0.07 ± 0.02. We also compared the results 
from O'Brian et al. (1991) with solar gf-values obtained by Thevenin (1989, 1990). We found that 
the mean difference in log(gf) is ±0.05 ± 0.02 when comparing 68 lines in common between these 
latest works. Similar offsets are found by Lambert et al. (1996). 

The gf-values for our Fe I lines were taken when possible from laboratory data in the following 
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order of priority: from O'Brian et al. (1991), from the Oxford Group, from May et al. (1974), and 
from Wolnik et al. (1971) corrected by —0.07 dex. In the cases where no laboratory data were 
available, we used solar gf-values from Thevenin (1989, 1990) corrected by +0.05 dex. 

Transition probabilities for the Fe II lines were taken from the solar analysis of Blackwell et al. 
(1980), Biemont et al. (1991), and from the semiempirical calculations of Kurucz (1993b). While 
restricting ourselves to lines present in our spectra, we compared the gf-values from Blackwell et 
al. (1980) and Kurucz (1993b) to the study of Biemont et al. (1991). There are 9 lines in common 
between Biemont et al. (1991) and Blackwell et al. (1980). The mean difference between the two in 
log(gf) is 0.14 ±0.02. A similar offset was found by Lambert et al. (1996). We found no significant 
difference between the results of Biemont et al. (1991) and Kurucz (1993b), since the 9 lines in 
common result in a mean difference of 0.03 ± 0.02. 

The gf-values for our Fe II lines were taken in the following order of priority: from Biemont 
et al. (1991), from Blackwell et al. (1980) corrected by 0.14 dex, and from Kurucz (1993b). 



Our sample contains many stars observed in mutiple orders, with many detectable absorption 
features in these high S/N spectra. For example, in the coolest M71 star in our sample, M71- 
1-45, 1407 absorption lines have been identified. A FORTRAN code to automatically search for 
absorption features and measure their equivalent width (W\), EWDET, was developed for this 
project. The code is available upon request to SVR. EWDET determines the continuum location 
of the HIRES spectra by fitting a curve to the spectra performing several iterations of point rejection 
above and below sigma levels given by the user. Then, EWDET identifies lines above the noise 
level defined by a factor of two of the continuum rms dispersion. Each of the identified lines is fit 
by a Gaussian profile and then the W\ are computed by the integration of the fitted Gaussian. The 
error in W\ is computed by adding quadratically the error at each point of integration, cij, times 
the step of the integration. The error at each point of integration is given by: 



where g(\i) is the Gaussian profile, P is the peak of the Gaussian, ap is the error in the peak of the 
Gaussian, A cen is the central wavelength of the Gaussian, <J\ cen is the error in the central wavelength, 
a is the dispersion of the Gaussian, a a is the error in the dispersion of the Gaussian, and a con t is 
the error in the continuum. The errors of the Gaussian parameters are from the covariance matrix 
of the Gaussian fit. The expression for the error in W\ is deduced by propagating the errors of the 
Gaussian parameters and assuming that the continuum level is equal to one (see Appendix). The 
fit by a Gaussian profile is reasonable even for the strongest lines we use, as shown in Figure 1, 
where the observed line is plotted with a solid curve and its corresponding Gaussian profile fit is 
shown with a dashed curve. 



2.2. Measurement of Equivalent Widths 
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The line list identified and measured by EWDET is then correlated to the line list with the 
atomic parameters to specifically identify the Fe I and Fe II lines. The detailed lists of W\ and 
gf-values will be given in the next paper. 

The spectral resolution, A/AA, of an echelle is fixed, unlike a low incidence, angle of low 
order grating spectrograph where AA is constant, independent of A. Hence a line of constant 
central depth, D, will have an equivalent width proportional to XD. We construct the relationship 
between XD and W\ for the Fe I lines to look for possible blends and for saturation effects. The 
XD versus W\ relation for three stars are plotted in Figure 2; M71-I is one of the most luminous 
and coolest stars of our sample, M71-G53476_4543 is a medium luminosity and temperature star, 
and M71-G53392_4624 is one of the faintest and hottest stars of the sample. We fit a second 
order polynomial to the XD—W\ relationship for each star, performing several iterations of point 
rejection above and below the two sigma level. The second order term is needed to account for line 
saturation at large W\ and the rejection of points is needed to eliminate blended lines. The second 
order fit is plotted as a solid line in Figure 2, and the points considered for such a fit are shown in 
black. 

For Fe I, we use two sets of lines. The first set, subsequently called "the weak line set" , contains 
those Fe I lines which are within two sigma levels of the XD—W\ fit, have W\ < 60 mA, and have 
errors less than a third of the W\. This set of lines produces a sample of the best weak Fe I lines 
with the most accurate W\ and the abundances derived using them will have a minimal dependence 
on the choice of microturbulent velocity. The second set, subsequently called "the good line set", 
consists of all the Fe I lines with errors less than a third of the W\ and with W\ computed from the 
fit for XD—W\ determined for each star. In future papers, the W\ of lines of other elements will 
be determined in the same manner as the Fe I lines of "the good line set" . This way a consistent 
comparison can be done among the resulting abundances without a restriction on the strength of 
the lines used. The "weak line set" is different for each star. Actually the weak lines for the stars 
near the main sequence are no longer weak for the cooler stars in our sample. We compare the 
results from 20 lines common to 15 stars over almost the whole range in T e ff with the results from 
the "weak line set"; there is no difference within the errors, nor a trend in T e ff. 

For Fe II, the W\ of the lines are also determined using the fit to the XD—W\ relation of the 
Fe I lines of "the good line set" . The Fe II lines follow the relationship determined from Fe I lines 
well, as shown in Figure 3. Additional Fe II lines, not picked up automatically, were measured by 
hand for the stars near the turn off and for M71-G53425_4612 and M71-G53457_4709. The set of 
W\ for the hand selected Fe II lines is computed from the XD—W\ relation of the Fe I lines, after 
determining their observed depth from the spectra. 

The number of Fe I lines, for both sets of lines, and the number of Fe II lines utilized in this 
analysis are listed in Table 1. 
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2.3. Spectroscopic Effective Temperature 

The effective temperature (T e ff) of a star can be determined spectroscopically by requiring the 
abundance to be independent of the lower excitation potential. This technique can be applied to 20 
of our stars where we have Fe I lines with enough range in lower excitation potential to do so. For the 
spectroscopic T e ff determination we are using "the weak line set" of Fe I lines, because its resulting 
abundance and spectroscopic T e ff will be only weakly dependent on the choice of microturbulent 
velocity. We find that the spectroscopic T e ff is in good agreement with the photometric T e ff 
derived in Paper I, as shown in Figure 4. The solid line in Figure 4 shows the ideal case when the 
spectroscopic and the photometric T e ff are equal. The scatter around the solid line is about 150 
K, which is comparable to the error of the photometric T e ff of 75 K for giants and of 150 K for 
dwarfs (Paper I), also shown in Figure 4. The scatter around the solid line is symmetric, not above 
or below, indicating the lack of systematic effects with the photometric temperatures. 

2.4. Microturbulent Velocity 

The microturbulent velocity (£) of a star can be determined spectroscopically by requiring the 
abundance to be independent of the strength of the lines measured as the equivalent width. We 
apply this technique for both sets of Fe I lines. The resulting £ and the [Fe/H] computed with it 
for "the weak line set" of Fe I lines are listed in Table 1. Only 18 of our stars have enough weak 
Fe I lines to derive £ spectroscopically. We also compute £ using all the good Fe I lines for 20 of 
our stars. The difference of the resulting [Fe/H] computed with the derived £ for the two set of 
lines is plotted with respect to T e ff in Figure 5. The mean difference is — 0.004± 0.011, hence, as 
expected, the [Fe/H] results from both sets of lines show a very good agreement. 

The relationship between £ determined for the set of all good Fe I lines and the photometric 
T e ff is shown in Figure 6. The solid line corresponds to a linear least squares fit of the data, 
excluding the red horizontal branch (RHB) stars, marked with circles. The best fit line is given by: 

£ = 3.30 -4.66 x 10" 4 xT eff 

The scatter around the solid line is about 0.2 km s _1 , which is a reasonable estimation of the error 
in £. For the rest of the analysis, we will use the set of all the good Fe I and Fe II lines, with £ 
computed from the £-T e ff fit. For the RHB stars we use a value of 1.61 km s _1 , which corresponds 
to the mean value determined for the three RHB stars. The microturbulent velocity used for our 
stars is listed in Table 1. 

3. RESULTS 

Given the stellar parameters from Table 1, we determined the iron abundance using the equiv- 
alent widths of all the good Fe I and Fe II lines identified in the HIRES spectra. We employ the 
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grid of stellar atmospheres from Kurucz (1993b) to compute the iron abundance using the four 
stellar atmosphere models with the closest T e ff an d log(g) to each star's parameters. The [Fe/H] 
listed in Table 1 is an interpolation of the results from the closest stellar model atmospheres to the 
appropriate T e ff and log(g) for each star. 

3.1. [Fe/H] from Fe I lines 

The results of [Fe/H] from Fe I lines are listed in column 9 of Table 1 and plotted against 
the photometric T e ff in Figure 7a. T e ff is used for the x-axis as a convenient parameter for 
characterizing the position of the stars in the color-magnitude diagram as it also ranks the stars 
in luminosity (except for the RHB stars). The errors listed in Table 1 correspond to the larger of 
the statistical uncertainty, given by the standard deviation of the iron abundance from different 
lines divided by the square root of the number of lines, or a minimum value of 0.03 dex. These 
errors are lower limits to the actual uncertainties in the abundances, since they do not include 
uncertainties due to the stellar parameters nor any systematic effects that might be present. We 
estimate the sensitivity of [Fe/H] derived from Fe I lines with respect to the stellar parameters in 
three cases 4250/1.0/1.0, 5000/2.5/1.0 and 5500/4.0/1.0, where the three numbers correspond to 
T e ff/log(g)/^. The results are listed in Table 2, where the range adopted for each parameter is 
representative of its uncertainty. Our determination of [Fe/H] from Fe I lines is most sensitive to 
errors in T e ff, which is less than ~ 0.1 dex for AT e jj of ± 100 K, and have a minimal sensitivity on 
the choice of metallicity of the model atmosphere grid for plausible changes in [Fe/H] (±0.2 dex). 

The solid line, shown in Figure 7a, is a linear fit weighted by the errors of [Fe/H] versus T e ff. 
The slope of the fit is (—0.8 ± 3.6) x 10~ 5 dex/K, which is consistent with [Fe/H] being constant, 
independent of T e ff (ie, of luminosity or equivalently position in the color- magnitude diagram). 
We divide our sample in four groups of stars: giant stars at or above the RHB, stars on the RHB, 
giant stars below the RHB, and main sequence stars near the turnoff. The mean [Fe/H] for each 
group is listed on Table 4. We found no significant difference in the mean [Fe/H] obtained from Fe 
I lines among the defined groups of stars. 

The mean [Fe/H] weighted by the errors of all 25 stars is —0.71 ±0.08, in very good agreement 
with earlier determinations (Cohen 1983; Gratton et al. 1986; Leep et al. 1987; Sneden et al. 1994). 

3.2. [Fe/H] from Fe II lines 

The determinations of [Fe/H] from Fe II lines are listed in column 11 of Table 1 and plotted 
against the photometric T e ff in Figure 7b. The errors listed in Table 1 corresponds to the statistical 
uncertainty or a value of 0.05 dex, whichever is larger. We estimate the sensitivity of [Fe/H] derived 
from Fe II ines with respect to the stellar parameters in the same manner as the sensitivity of [Fe/H] 
from Fe I lines. The results are listed in Table 3, where the range adopted for each parameter is 



-8- 



representative of its uncertainty. We see a stronger sensitivity on the stellar parameters from the 
Fe II lines than from the Fe I lines. The [Fe/H] determination from Fe II lines is most sensitive to 
the systematic error (note that the internal uncertainty in log(g) is < 0.1 dex) in log(g), as well 
as to T e ff among the coolest M71 giants. The sensitivity on the choice of metallicity of the model 
atmosphere grid is small for reasonable changes in metallicity. 

The solid line, shown in Figure 7b, is a linear fit weighted by the errors of [Fe/H] versus T e ff. 
The slope of the fit is (+3.1 ± 5.2) x 10~ 5 dex/K, which is consistent with [Fe/H] being constant, 
independent of T e ff. We found that there is no significant difference in the mean [Fe/H] obtained 
from Fe II among stars from different luminosity groups, listed in Table 4. The mean [Fe/H] 
weighted by the errors is —0.84 ± 0.12, in very good agreement with our result from Fe I lines and 
earlier determinations (Cohen 1983; Gratton et al. 1986; Leep et al. 1987; Sneden et al. 1994). 

3.3. NLTE effects 

The iron abundance could be affected by departures from LTE. The main NLTE effect in late- 
type stars is caused by overionization of electron donor metals by ultraviolet radiation (Auman & 
Woodrow 1975). Recently, Gratton et al. (1999) and Thevenin & Idiart (1999) studied NLTE effects 
in Fe abundances in metal-poor late-type stars. Gratton et al. (1999) found that NLTE corrections 
for Fe lines are very small in dwarfs of any T e ff, and only small corrections (< 0.1 dex) are expected 
for stars on the red giant branch. Thevenin & Idiart (1999) found that NLTE corrections become 
less important as [Fe/H] increases, being less than 0.1 dex for stars with [Fe/H]> —0.75 dex, and 
that ionized lines are not significantly affected by NLTE. 

One way to explore possible NLTE effects present in our data is by comparing the results 
from Fe I and Fe II lines. The difference between [Fe/H] from Fe II and Fe I lines is plotted in 
Figure 8 against T e ff. The solid line is a linear fit weighted by the errors. The slope of the fit is 
(+2.0 ± 8.2) x 10" 5 dex/K, which is nearly flat. The mean difference is —0.13 + 0.18. We conclude 
that NLTE effects are negligible in our iron abundance determination, as expected from results of 
earlier studies (Gratton et al. 1999; Thevenin & Idiart 1999). 

4. DISCUSSION 

Our [Fe/H] abundance results involve for the first time a wide luminosity sample of stars, 
which includes at the same time stars from luminous giants to stars near the turnoff. We find that 
the [Fe/H] abundance, from both Fe I and Fe II lines, is independent of T e ff, and equivalently 
luminosity. 

Our result is in agreement with the work of Gratton et al. (2001). They present abundances 
from high dispersion spectra from the VLT of stars in NGC 6397 and NGC 6752. They found that 
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the [Fe/H] obtained for stars at the base of the subgiant branch agrees within a few percent with 
the [Fe/H] obtained for stars at the main sequence turnoff, and further compare this value with 
analysis of the RGB stars in this cluster by other groups. Note that the luminosity range of the 
sample presented in our work is several orders of magnitude wider than the luminosity range of 
Gratton et al. (2001) 's sample. 

Our results, and those of Gratton et al. (2001), appear to be in disagreement with inhomo- 
geneities in [Fe/H] found earlier by King et al. (1998). They obtained [Fe/H]=— 2.52 dex for a 
sample of subgiant stars in M92, which is a factor of two smaller than [Fe/H] measurements using 
red giants in the same cluster (Cohen 1979; Sneden et al. 1991). King et al. (1998) compare their 
result for the M92 subgiants with analysis of RBG stars by other groups, who may have determined 
the stellar parameters and performed the abundance determinations in a different way. This pos- 
sible difference in the analysis of the giant and subgiant sample may account for the difference in 
[Fe/H] found by King et al. (1998) or perhaps the determination of the stellar parameters by King 
et al. (1998) is flawed. Our result, on the other hand, is robust, because we have determined both 
the stellar parameters and the Fe abundance in a homogeneous and consistent manner for all our 
stars. 

Gratton et al. (2001) also found that NGC 6397 is homogeneous in both O and Fe, while 
an O-Na anticorrelation is present among unevolved stars in NGC 6752, which is very difficult to 
explain by the deep mixing scenario. Lines from many additional species, including O, Na, Mg, Ti, 
Sc among others, are observed in our HIRES spectra. We plan to present additional information 
in the matter of light, iron-peak, and heavy elements in the near future. 

5. CONCLUSIONS 

We present results of a high dispersion analysis of Fe I and Fe II lines to obtain [Fe/H] for 
25 members of the Galactic globular cluster M71. Our sample of stars includes 19 giant stars, 3 
horizontal branch stars, and 3 stars near the turnoff. Our conclusions are summarized as follows: 

• The [Fe/H] obtained from Fe I lines agrees very well with the [Fe/H] obtained from Fe II 
lines. 

• The mean [Fe/H] obtained from Fe I and Fe II lines of all 25 stars is in good agreement with 
earlier determinations. 

• The [Fe/H] obtained from both Fe I and Fe II lines is independent of T e ff, and equivalently 
luminosity. 

• No difference is found among the mean [Fe/H] from giant stars located at or above the RHB, 
RHB stars, giant stars located below the RHB and stars near the main sequence turnoff. 
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In the near future, we will present the result from the analysis now underway of additional 
elements. 
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A. Errors in the Equivalent Widths 

The equivalent width of a line is defined as: 
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where, g(X) is the flux of the spectrum, cont(X) is the continuum level, and AA is the step that 
can be made arbitrarily small to better approximate the integral. In this case, g(X) is the Gaussian 
profile of the line, given by: 
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where P is the peak of the Gaussian, A cen is the central wavelength, and a is the dispersion of the 
Gaussian. If the equivalent width is computed through out the summation then the error in the 
equivalent width, a w will be the quadratic summation of the errors in the individual points times 
the step: 
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If we assume that the continuum level is equal to one then: 
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Since the Gaussian depends on three parameters: 
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Table 1. Stellar Parameters for the M71 Sample 



ID a T eff log(g) £ b Npei b [Fe/H]pei b £ c N FcI c [Fe/H] FcI c N Fen d [Fe/H] FeII d 

(K) (km/s) (km/s) 



1-45 


3950 


0.90 


1.48 


59 


-0.74 


± 


0.03 


1.46 


187 


-0.78 


± 


0.03 


6 


-0.61 


± 


0.11 


I 


4150 


1.00 


1.00 


67 


-0.69 


± 


0.03 


1.37 


186 


-0.76 


± 


0.03 


6 


-0.88 


± 


0.07 


1-66 


4250 


1.35 


1.80 


68 


-0.83 


± 


0.03 


1.32 


179 


-0.71 


± 


0.03 


6 


-0.86 


± 


0.09 


1-64 


4200 


1.35 


1.57 


61 


-0.76 


± 


0.03 


1.34 


187 


-0.74 


± 


0.03 


5 


-0.73 


± 


0.09 


1-56 


4525 


1.60 


0.81 


25 


-0.63 


± 


0.04 


1.19 


127 


-0.56 


± 


0.03 


2 


-0.83 


± 


0.14 


1-95 


4550 


1.65 


1.00 


79 


-0.68 


± 


0.03 


1.18 


184 


-0.67 


± 


0.03 


8 


-0.87 


± 


0.05 


1-81 


4550 


1.75 


1.50 


77 


-0.71 


± 


0.03 


1.18 


180 


-0.64 


± 


0.03 


6 


-1.05 


± 


0.05 


1-1 


4700 


2.05 


0.89 


55 


-0.62 


± 


0.03 


1.11 


134 


-0.67 


± 


0.03 


5 


-0.86 


± 


0.05 


l-80 c '^ 


5300 


2.45 




28 


-0.68 


± 


0.04S 


1.61 


71 


-0.69 


± 


0.03 


5 


-0.92 


± 


0.05 


l-87 e 


5300 


2.45 


1.62 


68 


-0.58 


± 


0.03 


1.61 


128 


-0.60 


± 


0.03 


9 


-0.83 


± 


0.05 


l-94 e 


5300 


2.45 




42 


-0.78 


± 


0.03S 


1.61 


94 


-0.77 


± 


0.03 


6 


-0.81 


± 


0.05 


1-60 


4900 


2.30 


0.70 


49 


-0.77 


± 


0.03 


1.02 


119 


-0.80 


± 


0.03 


6 


-0.70 


± 


0.05 


1-59 


4600 


2.30 


1.50 


62 


-0.79 


± 


0.03 


1.16 


141 


-0.78 


± 


0.03 


5 


-0.70 


± 


0.05 


G53476.4543 


4900 


2.65 


0.83 


94 


-0.64 


± 


0.03 


1.02 


174 


-0.68 


± 


0.03 


7 


-0.84 


± 


0.05 


2-160 


5100 


2.70 


1.10 


68 


-0.59 


± 


0.03 


0.92 


145 


-0.54 


± 


0.03 


5 


-0.97 


± 


0.08 


G53447.4707 


5175 


2.75 


1.35 


90 


-0.62 


± 


0.03 


0.89 


155 


-0.57 


± 


0.03 


7 


-0.86 


± 


0.05 


G53445.4647 


5050 


2.85 


0.54 


50 


-0.54 


± 


0.03 


0.95 


112 


-0.65 


± 


0.03 


6 


-0.85 


± 


0.05 


G53447.4703 


5000 


3.00 


0.90 


62 


-0.72 


± 


0.03 


0.97 


125 


-0.77 


± 


0.03 


4 


-0.80 


db 


0.05 


G53425.4612 


5150 


3.15 


1.40 


36 


-0.77 


± 


0.03 


0.90 


80 


-0.73 


± 


0.03 


2 h 


-0.91 


± 


0.08 


G53477.4539 


5150 


3.15 




56 


-0.66 


± 


0.03 g 


0.90 


119 


-0.70 


± 


0.03 


5 


-0.90 


± 


0.05 


G53457.4709 


5200 


3.35 


1.24 


58 


-0.78 


± 


0.03 


0.88 


93 


-0.78 


± 


0.03 


5 h 


-0.76 


± 


0.11 


G53391.4628 


5100 


3.35 




55 


-0.74 


± 


0.03S 


0.92 


106 


-0.84 


± 


0.03 


5 


-0.81 


± 


0.07 


G53417.4431 


5800 


4.05 




19 


-0.66 


± 


0.03S 


0.60 


38 


-0.68 


± 


0.04 


3 h 


-0.61 


± 


0.12 


G53392.4624 


5800 


4.05 




23 


-0.81 


± 


0.04S 


0.60 


36 


-0.81 


± 


0.03 


3 h 


-0.66 


± 


0.08 



Table 1 — Continued 



ID a 


T ef.f 


log(g) t 


N FeI b [Fe/H] FeI b 


e 


N Fe i c [Fe/H] Fe i c 


N Fe „ d [Fe/H] Fe „ d 




(K) 


(km/s) 




(km/s) 






G53414.4435 


5900 


4.15 


5 -0.82 ± 0.12« 


0.55 


13 -0.83 ± 0.04 


2 h -0.58 ± 0.17 



identifications are from Arp & Hartwick (1971) or are assigned based on the J2000 coordinates, rh rm rs.s dd dm dd 
becoming Grmrss_dmdd. 

b Set of weak Fe I lines. 

c Set of all good Fe I lines. 

d Set of all good Fe II lines. 
e RHB star. 

f Appears to show rotation (Paper I). 

g Computed with £ from the set of all good Fe I lines. 

h Includes additional Fe II lines selected by hand. 
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Table 2. Sensitivity of [Fe/H]p e i on Stellar Parameters 





+ 100 K 


Alog(<?) 
+ 0.2 dex 


+ 0.2 km s" 1 


A[Fe/H] 
+ 0.2 dex 


4250/1. 0/1. a 


+0.04 


+0.02 


-0.08 b 


-0.03 


5000/2.5/1.0 a 


+0.09 


+0.01 


-0.06 b 


-0.01 


5500/4.0/1.0 a 


+0.08 


+0.02 


-0.03 b 


-0.01 



a T e// /log( ff )/e 

b This is for the set of good Fe I lines. It is smaller by a factor of 3 
for the set of weak Fe I lines. 



Table 3. Sensitivity of [Fe/Hjpeii on Stellar Parameters 





^T ef f 
+ 100 K 


AlogG?) 
+ 0.2 dex 


+ 0.2 km s" 1 


A[Fe/H] 
+ 0.2 dex 


4250/1.0/1.0* 


-0.12 


+0.11 


-0.04 


-0.07 


5000/2.5/1.0 a 


-0.02 


+0.09 


-0.03 


-0.04 


5500/4.0/1.0 a 


-0.03 


+0.08 


-0.02 


-0.03 



a T e// /log(<7)/£ 



-17- 



Table 4. [Fe/H] for Each Group of Stars 



Star group 


< V > 


N s tars 


< [Fe/H] FeI > 


< [Fe/H] FeII > 


RGB at or above RHB 


13.46 


10 


-0.71±0.07 


-0.83±0.12 


RHB 


14.50 


3 


-0.68±0.07 


-0.86±0.05 


RGH below RHB 


15.92 


9 


-0.69±0.09 


-0.85±0.06 


MS TO 


17.76 


3 


-0.78±0.06 


-0.64±0.13 
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Fig. 1. — Strongest observed Fe I lines for M71-1-45. The observed lines are plotted with a solid 
line, and the corresponding Gaussian profile is plotted with a dashed line. The W\ of each line is 
indicated in the botton left corner of each pannel. 
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Fig. 2. — Central depth times wavelength versus W\ relations for M71-I (one of the most luminous 
and coolest stars in our sample), M71-G53476_4543 (a star of medium luminosity and temperature), 
and M71-G53392_4624 (one of the faintest and hottest stars in our sample). The solid curve is a 
second order fit obtained after several iterations of rejection of points deviating by 2a or more. The 
points used in the fit are shown in black. 
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Fig. 3. — Central depth times wavelength versus W\ relations for M71-I (one of the most luminous 
and coolest stars in our sample), M71-G53476_4543 (a star of medium luminosity and temperature), 
and M71-G53392_4624 (one of the faintest and hottest stars in our sample). Solid squares denote 
identified Fe II lines and the solid curve is the second order fit obtained for Fe I lines and shown in 
Figure 2. 
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Fig. 4. — Photometric T e jj versus spectroscopic T e ff for the M71 sample. The solid line indicates 
the ideal case when the photometric and spectroscopic T e ff have the same value. The scatter around 
the solid line is about 150 K. The only RHB star for which T e ff can be determined spectroscopically 
is marked with an open circle. 
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Fig. 5. — The difference between [Fe/H] computed with the derived £ for the set of weak Fe I lines 
and with £ from the set of all good Fe I lines is plotted with respect to T e f f . The solid line indicates 
equality. The only RHB star for which £ can be determined spectroscopically is marked with an 
open circle. 
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Fig. 6. — £ determined for the set of all good Fe I lines is shown as a function of T e ff. The solid line 
is the linear fit weighted by the errors. The RHB stars are excluded from the fit and are marked 
with circles. The scatter around the solid line is about 0.2 km s _1 . 
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Fig. 7. — [Fe/H] from Fe I (upper panel) and Fe II (lower panel) against photometric T e ff. The 
solid lines are linear fits weighted by the errors. In both cases, [Fe/H] shows no dependence with 
T e ff. The dashed lines indicate the mean [Fe/H] with their respective error plotted as an error bar 
at 3925 K. Note that < [Fe/H] FeI >= -0.71 ± 0.08 and < [Fe/H] FeII >= -0.84 ± 0.12. The 
RHB stars are marked with a open circles. 
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Fig. 8. — Difference between [Fe/H] from Fe I and Fe II against T e ff- The solid line, which is nearly 
flat, is a linear fit weighted by the errors. The dashed line indicates the mean difference with its 
respective error plotted as an error bar at 3925 K. Note that the mean difference is —0.13 ± 0.18. 
The RHB stars are marked with a open circles. 



